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Nomenclature 


A 

Cut-off frequency in the spatial frequency domain 

c 

Path of radiation 

CAT 

Computer Aided Tomography 

CBP 

Convolution backprojection 


Inherent error 

tD 

Discretization error 

ERL 

Average error as defined in 3.2.2 

ERMIN 

Minimum error of the reconstructed image 

ERMAX 

Maximum error of the reconstructed image 

F 

Spatial frequency 

FBG 

Fan beam geometry 

NRAY 

Number of rays in a view 

NVIEW 

Number of views 

p{s; 9) 

Projected data 

PBG 

Parallel beam geometry 

R. 

Cui^off frequency 

W{F) 

Window or filter function 

As 

Distance between two consecutive parallel lines 

4>) 

Point dependent attenuation coefficient of the material 

p 

Density of the given material 

9 

Angle between the Y-axis and the given ray in PBG 
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Abstract 


CT is being used in various areaa. Of the various algorithms reported so fax, e.g. 
ConvdutioQ BackprojectiocL (CBP), ART, Direct Fourier Inversion etc. the CBP 
algorithm has gained wide acceptance and popularity owing to its efficiency and ease 
of application. Some error estimateB for the CBP algorithm have been reported 
by Natterar and MunsM et al, which incorporate the assumption of the projection 
data being essentially band limited. Due to finite frequency cut-off incorporated in 
implementation of CBP algorithm, the approximation obtained is having inherent 
error, A practical implementation of CBP method also introduce, apart from 
eji, the discretization error arising due to discrete implementation of the convolution 
and backprojection integral. A numerical study of these errors for recoustruction of 
simulated images, using CBP was carried out. The results show that they axe in tune 
with the theoretical predictions. Errors are found to be of the order 0(1/ Re). 
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Chapter 1 
Introduction 


The development of imaging systems based on the concept of computerized tomogra- 
phy has revolutionised the area of medical imaging. The potential of CT scanners is 
adso being harnessed in the field of non destructive evaluation, at a large scale. This 
technique has proved its worth not only in the detection of cancerous tissues, but also 
in the area of void fraction measurment of two phase flow. 

Problem of determination of cross sectional distribution of any property arises in 
many fields. Direct methods to solve the problem are not possible m many of the 
cases like non destructive testing, two phase flow, and medical imaging etc. Indirect 
methods axe employed in such cases. 

Projection data (refer chapter 2) can be obtained without much problem, by 
making use of x-ray, gama ray, or ultrasonic source detector system. Now the major 
problem to be dealt with is reconstruction of image from the projection data. Radon 
[1] paved the way to solve this hurdle in 1917. He proved that any arbitreiry function 
could be recovered from its set of line integral taken along various chords and various 
directions, using Radon inversion formula. But this could not be implemented due to 
inherent mathematical complexities. BraceweU [2] reported an application of CT in 
radio astronomy and Cormack[3] derived an inversion formula in 1963. These two were 
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more closer to implementation tkan the Ration inversio formula™ Bracewell ^ Riddle 
[4] made use of convolutions in 1967, which accelerated the CT computations. The 
use of RL filter by Ramachandran &: Lakshminaiayan[5] further accelerated the CT 
methodology. But it was Hounsfield[6] who pioneered the way to practical application 
of CT in 1973. 

Some work hzis been done at IIT Kanpur, which deals with many fundamental 
problems in the field of CT[7]. This involves void firaction measurement using algo- 
rithms of gamma-ray tomography[8], estimation of error for tomographic inver8ion8[9], 
and developing new filters for image reconstruction[10,ll]. 

This study involves reconstruction of given images, using CBP algorithm, for 
different number of rays and studying the error in the reconstructed images. 



Chapter 2 


Preliminaries 


When a beam of photon is incident on matter, it can interact with matter in many 
ways. Some of the photons are scattered, some undergo absorption, and rest are trans- 
mitted. The transmitted beam is said to be attenuated. The attenuation depends on 
thickness of the object and energy of the incident beam. 

Single beam mono-energetic radiation attenuation phenomenon in a plane can be 
represented as follows : 

N = No exp[- ^{r, 4>) d/] (2.1) 

where 

N is no of incident photons, 

No is no of detected photons, 

IJ is Hnear attenuation coefficient for given energy, 

I is the integration variable, 

c chord, along which function is to be integrated, and 
r, <i> axe cylindrical co-ordinates. 

The value of n for a gjven material depends on the type of radiation and energy of . 
the radiation. In the equation (2.1) energy dependence of /i is ignored, n is considered 
to be a function of r and 4> only, as the radiation is restricted to a plane only. 
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The equation can be rewritten as : 

P = ( 2 . 2 ) 

where 

p = ln(JV„/iV), (2.3) 

Now we have to find the value of ^^(r, (f) for getting the density distribution. 
Radon[l] showed that n cam be recovered from a set of p-values measured eilong 
various strips, using following formula : 

= -— lirn / - mi(r coe <f> + Tsin<j> + q, <j>) d4> dq^ (2.4) 

e—O Jt q JO 

where mj denotes partial derivative of /i(r, <j}) with respect to r. Rigorous proof of 
this formula is discussed by Herman[12]. The above written Radon formxila has many 
implementation problems. In CT we have only a finite set of measurements, which 
are not sufficient to reconstruct the image accurately. Moreover the Radon formula 
is very sensitive to some inaccuracies which are inherent in CT like beam hardening, 
photon statistics, detector inaccurauries etc. 


2.1 Data Collection Modes 

Projection data are required for CT, collected by an array of radiation detectors for 
reconstruction of the function /i(r, 0). Various data collection modes axe discussed 
by Herman[12]. Two popular modes of data collection are termed as paiallel-beeun 
geometry (PBG) and fan-beam geometry (FBG). 

2.1.1 Parallel-Beam Geometry 

PBG mode CT scanners have a no. of uniformly spaced source-detector pairs. The 
object is mounted at the centre, around which axe located the source-detector pairs. 



CHAPTER 2. PRELIMINARIES 


5 


The object can be rotated to ^ve different values of angle d. For eaxi. value of 0 a set 
of projection data is collected. Figure (2.1) depicts PBG configuration. 

The line SD represents the path of the data ray. The perpendicular distance of 
the ray from the origin is denoted by s. Several SD pairs collect the data p for a given 
9, This set of p is known as a ‘projection’, which are collected for different views ( at 
different 6 ). This data is denoted by p(s; 6). 

2.1.2 Fan— Beam Geometry 

In this mode, a single source is viewed by several detectors, simultaneously. The 
source angle is denoted by a and the detector angle by The readings are taken 
at different values of to get the data, denoted by g{c,^) or h{l,P). Here, I is 
the perpendicular distance from the origin to the particular ray. This configuration, 
depicted in Fig-2.2, is widely used in tomography for medical purposes.The FBG 
mode can be implemented for a rotating object system as well. The data collection 
methodology remains the same as that for the PBG case. 


2.2 Mathematical Formulae for Tomographic In- 
version 

The mathematical methods used to retrieve a two dimensional function from its pro- 
jection data, known as Tomographic inversion methods, are generally categorized as 
the ‘transform methods’ and the ‘algebraic reconstruction techniques’- The transform 
methods are based on mathematical formulae while the latter are iterative compu- 
tational techniques. Algebraic reconstruction is useful when the no. of projection 
data available is small. Whereas transform methods give accurate results only with 
sufficient projection data. In the following sections, ‘transform” methods are briefly 
summarised for the two data collection geometries. 
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2.2.1 Parallel Beam Formulation 

The reconstniction process of a two-dimensional function from its projections (line 
integrals ) involves Fourier transforms for parallel beams and Hilbert transforms for 
divergent beam data. A very convenient transformation from FBG to PEG avoids the 
use of the Hilbert transform. The ‘central-slice’ theorem (also called the projection 
theorem ) states that the one dimensional Fourier transform of the projection data^ 
p(s; 9), with respect to the first variable s, is equed to the two-dimensional Fourier 
transform of the object function /. Mathematically, 

/(f;9)=p(f;9) (2.5) 

Taking the inverse Fourier transform of (2.5) we get 

/(r, ^)= / f p{F]d) exp[—i27cFrcos{d — <f}]\F\ dF dd (2.6) 

Jo CO 

For a given 9, 


p(F;9)= f p(s;^)exp[— :27rFs]ds. (2.7) 

J-R 

Elquation (2.6) requires continuous projection data for all values of s and 9. For 
computational feasibility, a filter function is incorporated in the Fourier domain, to 
enable a finite cutoff frequency. This necessity ^uises because the limits on R vary 
from — oo to -foo, which introduces divergence. This filter function, W{F), vanishes 
for I J"! greater thain A, the cut-off frequency. Equation (2.6) is transformed to 

f{r,<f>)= / p{F]9) exp[-i2TFrcos{9^<f>)]W{F)\F\dFd9. (2.8) 

Jo J— CO 

The reconstruction is approximate as all the higher frequencies have been elimi- 
nated. Besides, as per the sampling theorem, the cut-off frequency and the sampling 
interval in the spatial domain are related by 


A > [l/(2As)]. 


( 2 . 9 ) 
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where As is the spacing between the rays. 

The band-limiting filter, introduced in electron micrography by Raunachandran & 
Lakshininarayan[5] is given by 


W{F) = 


1, \F\<A, 
0, \F\>A. 


The sine filter used by Shepp & Logan[13] is given by, 


( 2 . 10 ) 


W{F) = 


[sm{'KF/2A)]/{tF/2A), jFj < A, 

0, IFl > A 

The filter used in this thesis is the Ramchandran filter. 


( 2 . 11 ) 


2.2.2 Fan— Beam Formulation 

The inversion formula for the FBG case was first derived by Herman and Naparstek[14] 
and is given by, 

“ 4^ / ~ 

where, 

D = distance of source firom reference origin, 

<x = angle of the data-ray with the reference ray, 

B = angle of the extreme data rays, ^ — source position, 
g{a, 0) = data for the ray represented by (<x, j0), 

D, = {ilu){d/d<T - d/d0), 

cd = tan~^[(rco8(;0 — if>))/{D + r^{0 — ^))], 

U = [(r co6(;5 - <f>)y + {D + rM0 " 
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amd the remainiiig vaxiables are same as for the PBG case (see Fig.(2.l). In the 
above equations U is the distance of the radiation source from (r, 4>), the point being 
reconstructed ,and o' is the angular displacement of the pairticulax data ray passing 
through that point (r, ^). A major difference ( compared to Eq.( 2.7), which derives 
the reconstruction formula for parallel beam geometry) is the computation of partial 
derivatives of the data, g. 

2.3 Convolution Back-Projection Algorithm 

Reconstruction algorithms based on the inversion formulae are ceiiled ‘transform 
methods’[15]. They are slow and require accurate interpolating schemes in computing 
the two-dimensional inverse Fourier transform. The introduction of convolution[4,5] 
eliminated the use of the Fourier transform and its inversion. For PBG , Eq.(2.8) can 
be written as 

/(r, (f>)= f I p(s; $) q(s' - s) ds d$, (2.13) 

•Jo 

where 

q{s) = W{F) IF’! exp(t’27rF’s) dF, (2.14) 

J — A 

and s' = rcos(d — <j>). 

Here q is known as the convolving function and is the inverse Fourier transform 
of the function, IV(F')|jF|, where W{F) is the filter function. The inner integral of 
Eq.2.12) is a convolution and the outer integral is called the back-projection. The 
CBP algorithm is the most widely used method of reconstruction employed by com- 
mercial CT scanners used in the area of medical imaging. 


Chapter 3 


Program Implementation and 
Description 


The CBP algorithm was used for reconstruction of 10 im^es and the error was 
calculated for each image for different no. of rays. The band limiting filter used for 
the image reconstruction was Ramchandran filter. 


3.1 Computer Implementation of the CBP Algo- 
rithm 

The problem of reconstruction from projections is to find f{x, y) from the gjven value 
of p{s‘, t). In practice, the problem may be stated as ; 

Given discrete projection data in the form of estimates of p for a finite number of 
rays, find a 2-D distribution, which is a reconstructed estimate of the unknown object 
In the case when p is seunpled uniformly in both s and 6, for N angles Ad apafft, 
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with each view having M equispaced rays As apart, we clefiiie[l5] 

Af+ = (M - l)/2 
3f- = -(M-l)/2 


Modd 


M even 


M+ = {M/2) - 1 
M- = -M/2 

In order to ensure that the collection of rays specified by 


|(rnAs, nA(9): < M~, 1 ^ n < 


covers the unit circle, we have, 


(3.1) 


(3.2) 


Ae = 7r/iV As = 1/Af+. 


(3.3) 


A reconstruction algorithm which can be implemented on. the digital computer 
is required to ev-aluate fjg{kAx, I Ay), which is a band limited approximation of the 
function to be reconstructed. Here K~ < ^ < K'^ and L~ < I < where k and I 
are the positions of the co-ordinates of the image pixel. The definition of their upper 
and lower limits is similcur to that of ‘m’. Thus the projected data from N views and 
M rays is to be used to construct an image oi K x L pixels. The back-projection 
integral is evaluated as follows : 

N 

fs{kAXf I Ay) o; Ad p{fs{kAx cos -b lAy sin 0„) (3.4) 

n=l 

For each angle the convolved values of p(s', 0„) for the KxL values of s', we can 
either have a separate convolution for every s' with the actual value of q{s' — mAs) 
at that point or we can evaluate p(mAs, 0^) only within the specified limits of m 
and then use interpolation. The latter approach is much faster and cheaper. These 
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operations axe represented by 


M+ 

Pc{m'As, $„) ~ As ^ p{TnAs, d„) q {{m' - m)) , 

M~ < m < M'^ 

— As^p(mAs,l9„) /(s' -- mAs) 


( 3 . 5 ) 

( 3 . 6 ) 


where I{s) is an interpolating function. A linear interpolating function, say /x(s), 
corresponding to linear interpolation between adjacent samples is 


h{s) = { 


^(l-ls|/As), l3|<As 
0, jj] < As 


(3.7) 


The other popular scheme is nearest neighbour interpolation. Both of them give 
satisfactory results. The choice depend on the information required by the usercite- 
herman. 


3.2 Errors Involved in Tomographic Inversion 

Implementation of CBP algorithm involves various errors, which can be categorized as 
inherent error and discretization error. The error estimation for tomographic inversion 
is discussed in detail by Munshi et al.[16]. 

3.2.1 Inherent error 

The inherent error is given by 

«/l = f{r, 4>) - /(r, 4>) (3.8) 

where /(r, <f>) and /(r, (f) are pven by equation 2.6 and 2.8. /(r, 4>) represents ap- 
proximation of /(r, <j>y This error arises due to finite cut-off of the Fourier frequency. 
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This error is precisely zero if the function /(r, <f>) happens to have a band limited 
projection data, and cut-off frequency is chosen such that it is at least equal to the 
Icvrgest frequency contained in /, and w{R) = 1, in the limi t [—Rc, R^. 

3.2.2 Discretization error 

The data obtained for reconstruction Me not continuous, but discrete. This causes 
error in the implementation of algorithm due to discretization of convolution and 
back projection integred. 

The total error in the reconstruction, at any point can be ^ven as sum of the 
two aforesaid errors, at that point. 


Sr ~ d" 


(3.9) 


The errors involved m the reconstruction of the simulated images are represented 
in the foUawing manner : 


Here, 


ERMIN 

ERMAX 

ERLl 

ERL2 

U 


nun(/-/) 

max(/-/) 

h'LElih - h)\ 


MP 

1 ' 
MP 


* } 


cd\ 


_ • J 


1/2 


NP 




» 3 


(3.10) 


MP represents the total number of pixels considered for evaluation of ERLl and 
ERL2, 


NP represents the total number of pixels considered for evaluation of E 



CHAPTER 3. PROGRAM IMPLEMENTATION AND DESCRIPTION 


15 


and, superscript CD impKes the value of the function after the 32 gray level coding. 


The errors, ERLl and ERL2 represent the usuzd /i and I 2 errors. The error 
E represents an index of visual quality of the reconstructed image. For a given 
window and fixed cutoff frequency R^, ERMIN and ERMAX give an indication of 
smoothness of the function concerned. 

The normalized errors corresponding to ERLl, ERL2, and E are defined as fol- 
lows: 


ERLIN 

ERL2N 

EN 


(ERLl) (MP) 
(ERL2) (MP) 


1, 




Zi E, /.5 


rCD 


(3.11) 


Here, / denotes the function value averaged over the entire region of interest. 



Chapter 4 


Results and Discussion 


This chapter deals with the data i.e. images used for reconstruction and the errors 
obtained on implementation of the CBP algorithm. 


4.1 Data Used 

The simulated images used for reconstruction are listed in Table 4.1. These images 
use a 32 gray level scheme to denote the density at different points. The ori^al 
images are shown in Figs. 4.1- 4.3. 


4.2 Results 

Reconstructed images for different number of rays are shown in the appendix. The 
various errors i.e. average errors li, ERMAX, and ERMIN are displayed in 
the Tables 4.2 - 4.5 for different number of rays. The plots depicting variation of 
ERMAX and ERMIN with Fourier frequency are shown in Figs. 4.4 - 4.23 and in 
the Appendix— B. The reconstructions with varying Rc indicate that error reduces in 
a linear fashion. It can be seen that order of ex ** 0{1/Rc). 
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4.3 Discussion 

As it is expected, the total error seems to be decreasing with the increase in number 
of rays. The total error consists of inherent error and discretization error. Both 
the errors 2 ire coupled, as the Fourier cut-off frequency is related to the number of 
rays. The nature of the curves obt^ed between error and Fourier frequency gives a 
fair reflection that the error decreases at a faster pzu:e initially, and after a Fourier 
frequency value 64, the error decre 2 ises at a slower rate. Although, quantitatively it 
is difficult to make any conclusive renwk, yet qualitatively it can be very well said 
that the error curves have two distinct regions. The first region has drooping kind of 
characteristic. Whereas the second region shows a linear decrease in the error values. 
The least squares fit for the higher Fourier frequency versus error is also shown in the 
plots (Figs. 4.4 - 4.23). This characteristic seems to be in tune with the theoretical 
predictions reported by Natterer[l7] and Munahi[9]. 
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Table 4,1: List of images used for reconstruction 


Si No. 

Image 

Image Key 

1 

SAT 

V 

2 

LIN 

W 

3 

FIN 

X 

4 

GIN 

Y 

5 

JET 

A 

6 

CRO 

B 

7 

THO 

C 

8 

TST 

D 

9 

BRN 

E 

10 

PIC 

Z 

11 

IND 

N 













Figure 4.2: Original simulation of TST, SAT, LIN and FIN 
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Tiible 4.2: Suinmary of errors in recoastructioa of images SAT, LIN and FIN 


Image 

Name 

Image 

Code 

NRAY 

ERMIN 

ERMAX 

ERLi 

ERL 2 

SAT 

V-1 

64 

-5.3244 

4.8050 

1.1386 

1.4398 


V-2 

126 

-2.1813 

1.9667 

0.3678 

0.4875 


V-^ 

256 

-1.0952 

1.0670 

0.1632 

0.2303 


V-4 

512 

-0.5420 

0.5735 

0.1106 

0.1384 


V-5 

1024 

-0.1210 

0.4017 

0.1563 

0.1631 

LIN 

W-1 

64 

-5.0369 

4.3424 

0.8991 

1.1740 


W-2 

128 

-2.0879 

2.1520 

0.3137 

0.4264 


W-3 

256 

-L0424 

1.1817 

0.1362 

1 

0.1932 


W-4 

512 

-0.4148 

0.6560 

1 

0.0889 1 

0.1141 


W-5 

1024 

-0.1308 

0.3845 

0.1210 

0.1309 

FIN 

X-1 

64 

-7.3287 

9.7300 1 

1.6054 

2.0348 


X-2 

128 

-3.1491 

2.6086 

i 

0.6067 

0.7524 


X-3 

256 

-1.3809 

1.2625 

0.2784 

0.3435 


X-4 

512 

-0.4684 

0.6519 

0.1439 

0.1791 


X-5 

1024 

-0.1489 

0-4301 

0.1383 

0.1591 
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Table 4.3: Summary of errors in reconstruction of images GIN, JET atnd CRO 



Image 

Code 

NRAY 

ERMIN 

ERMAX 

ERLi 

ERIn 

GIN 

Y-1 

64 

-11.8286 

12.6952 

1.9095 

2.5491 


Y-2 

128 

-05.2623 

03.3543 

0.6273 

0.8680 


Y-3 

256 

-01.8502 

01.8295 

0.2644 

0.3980 


Y-^ 

512 

-00.7914 

00.8206 

0.1300 

0.1933 


Y-5 

1024 

-00.3833 

00.4918 

0.1101 

0.1408 

JET 

A-1 

64 

-9.1101 

7.4484 

1.0398 

1.3839 


A-2 

128 

-2.8679 

3.3377 

0.3538 

0.4990 


A-3 

264 

-0.8881 

1.4410 

0.1522 

0.2202 


A-4 

512 

-0.3940 

0.8061 

0.0876 

0.1205 


A-5 

1024 

-0.1078 

0.5027 

0.1048 

0.1205 

CRO 


64 

-5.6963 

4.4814 

0.7492 

i 

0.1008 



128 

-2.4807 

2.2000 

0.2720 

0.3913 


B-3 

256 

-1.1586 

0.7974 

0.1214 

0.1762 


B-4 

512 

-0.5655 

0.4923 

0.0614 

0.0874 


B-5 

1024 

-0.2486 

0.2294 

0.0469 

0.0581 
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Table 4.4: Summary of errors in reconstruction of images THO, TST and BRN 


Image 

Name 

Image 

Code 

NRAY 

ERMIN 

ERMAX 

ERL^ 

ERL 2 

THO 

C-1 

64 

-7.1519 

7.2522 

1.3028 

1.7130 


C-2 

128 

-2.9083 

2.4679 

0.4409 

0.6041 


C-3 

256 

-1.3023 

1.4420 

0.1951 

0.2768 


C-4 

512 

-0.8100 

0.6728 

0.1047 

0.1441 


C-5 

1024 

-0.3289 

0,4173 

0.1012 

0.1153 

TST 

D-1 

64 

-6.7023 

5.1341 

0.7994 

1.0601 


D-2 

128 

-3.6107 

2.6479 

0.2243 

0.3767 


D-3 

256 

-1.7692 

1,1559 

0.0912 

0.1752 


D-i 

512 

-0.7808 

0.5594 

0.0437 

0.0836 


D~5 

1024 

-0.3132 

0.3327 

0.0368 

0.0523 

BRN 

E-l 

64 

-8.4734 1 

9.7270 

1.4381 

1.9822 


E-2 

128 

-2.7575 : 

5.1121 1 

0.5012 

0.7097 


E“3 ! 

264 

-2.2170 

1.5601 

0.2367 

0.3461 


E-4 

512 

-0.7367 

0.8860 

0.1267 

0.1733 


E-5 

1024 

-0.2460 

0.4655 

0.1282 

0.1437 
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Table 4.5: Sunmiary of errors in reconstrnction of images PIC and IND 


Image 

Name 

Image 

Code 

NRAY 

ERMIN 

ERMAX 

ERLi 

ERL 2 

PIC 

Z-1 

64 

-9.7068 

9.7919 

1.0897 

1.7531 


Z-2 

128 

-2.6595 

2.6963 

0.3678 

0.6124 


Z-3 

264 

-1.1953 

1.2619 

0.1581 

0.2841 


Z-4 

512 

-0.8794 

0.7108 

0.0729 

0.1346 


Z-5 

1024 

-0.3026 

0.3683 

0.0489 

0.0747 

IND 

N-1 

64 

-8.7730 

8.3720 

1.5107 

1.9961 


N-2 

128 

-3.6732 

3.9831 

0.4918 

0.7304 


N-3 

264 

-1.8676 : 

j 

1.9215 

0.2311 1 

0.3411 


N-4 

512 

-0.7342 

1-1708 

0.1129 

0.1706 


N-5 I 

1024 

-0.3726 

0.6003 

0.1019 

0.1227 
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Figuie 4.4: Maximum error variation at higher frequencies for SAT 
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Figure 4.7: Minimum error variation at lugher frequencies for LIN 




CHAPTER 4. RESULTS AND 


fi 
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Figure 4.10: Maximum error variation at higher frequencies for GIN 
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Figure 4.11: Minimum error variation at higher ftequencies for GIN 
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Figuie 4.14: Maximum error variation at higher frequencies for CRO 
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Figure 4.18: Maximum error vaadatiou at higher irequencies for TST 
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Figure 4.22: Maximum error variation at higher frequencies for PIC 
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Conclusions and 
Recommendations 


N 

The application of Computer Tomography in a variety of fields ranging from medical 
imaging to astronomy demands a careful analysis of various errors involved in recon- 
struction of the projection data as different applications demand different precision 
and accuracy. The various errors involved in the reconstruction of objects are : error 
due to beam hardening, photon statistics, focal spot of the source and finite detector 
width (they are also referred as Partial Volume effect), organ motion effect etc. And 
of course, the discretization and inherent error. 

Beauty of using computer simulation lies in the fact that effect of these errors 
can be studied separately and they can be coupled to make a judicious choice of 
various parameters affecting the error of reconstruction. The importance of this 
investigation is obvious because this not only helps to decide number of rays to be 
used £31 particular appKcation, but also depicts the limitaticms of the algorithm nsecL 
Hence this investigation could be of immense przu:tical application for deciding various 
parameters in design and development of CT scanners for vanons purposes 

It is clear that the error value is reasonably small with number of views 1600 and 


46 



CHAPTER 5. CON CL USIONS AND RECOMMENDATIONS 


47 


mimber of rays 256. 

Following recommendations can be made regarding fuitber investigations: 

1. Discretization and inherent errors are coupled here and hence a scope to devise 
a methodology to separate the two and analyse their effect separately is there, 

2. The studies can be done by coupling the internal errors(such as discretization 
and inherent errors ) and external errors ( such as photon statistics and Beam 
hardening etc.) for particular applications. This can be of immense practical 
application. 

3. Effect of various filters for optimal mimber of rays and views on the error can 
be studied. 
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Figure A.2: Reconstructed Images 
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Figure A.5: Reconstructed Images 
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Figure A. 6: Reconstructed Images 
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Figure A.8: Reconstructed Images 
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Figure A. 12: Reconstructed Images 


Appendix B 


Error Variation with Frequency 
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Figure B.l: Error Variation With Fourier Frequency for SAT 
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PigTire B*2: Error Variation Witt Fourier Frequency for SAT 
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Figure B-.4: Error Variation With Fourier Frequency for LIN 
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Figure B*6: Error Vanaiiioii Witt Poimer Frequency for FIN 
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Figure B.7: Error Variation With Fourier Frequency for GIN 
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Figure B.9: Error Variation With Fourier Frequency for JET 
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Figure B. 13 : Error Variation Witt Fourier Frequency for THO 
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Fignre B.15: Error Variation Witk Fourier Frequency for TST 




78 


10.0 n 



FOURIER FREQUENCY 


Fignxe B.17: Eixrox Vaxiatioii Wiik Fourier Frequency fcr BH-N 


<>j - 




Wiafcion With Fourier Frequency for BRN 






Bibliography 


[1] J. Radon, Uber die Bestimmung von Funktionen durch ihre Integreilwerte langs 
gewifiser Maanigfaltigkeiten, Berichie Sachsische Akademie der Wissenschafien 
Leipzig, Matk.-Phys. Kl. 69, (1917) pp.262-267. 

[2] R. N. Bracewell. Strip integration in radio astronomy, Aust. J, Phys., 9, (1956) 
pp. 198-217 

[3] A. M. Cormack. Representation of a function by its line integrals with some 
radiological applications. J. Appl. Phys., 35, (1964) pp. 195-207 

[4] R. N, Bracewell, and A. C. Riddle. Inversion of fan beam scans in radio astron- 
omy. Astrophy. J., 150, (1967) pp.427-434 

[5] G. N. Ramachandrcin, and A. V. Lakshminarayan. 3-D reconstruction from radio- 
graphs and electron micrographs : Application of convolution instead of Fourier 
transforms. Froc. National. Academy. Science , USA, 68, (1971) pp.2236-2240 

[6] G. N. Hounsfield. Computerized transverse axial scanning tomography. Part I : 
Description of the system- Br. J. Radiol., 48, (1973) pp.1016-1022 

[7] P. Munshi, A review of computerized tomography with application to two-phase 
flows, Sadhana, 15, (1990) pp.43-55. 


81 





BIBLIOGRAPHY 


22 


[8] P. Munshi, Two-ph ase Flow Studies in the Bubbly Flow Regime Usiag a Scanning 
Gamma-Ray Densitometer, Master of Science Thesis, Ohio State University, 
(1979). 

[9] P. Munshi, Error estimates for the convolutioii backprojection algorithm in com- 
puterized tomography, Doctor of Philosophy Thesis, IIT Kanpur (1988). 

[10] R. K. S. Rathore, P. Munshi, P. Arora, S. D, Malik, A. K. Vaish, K. S. Singh, 
and U. Singh, A new non-Fourier tomographic filter for image reconstruction, 
Nucl. Technoi, 85, (1989) pp.346-349 

[11] P. Arora, P. Munshi, and R, K. S. Rathore, Higher order tomographic filters for 
non destructive purposes, Nucl. TecknoL, 83, (1988) pp.228-230. 

[12] G. T. Herman, Image Reconstruction from Projections : The Fundamentals of 
Computerized Tomography, Academic Press, New York (1980) 

[13] L. A. Shepp and B. F. Logan, The Fourier reconstruction of Head Section, IEEE 
Trans. : Nucl Sci., bf NS-21, (1974) pp.21^3. 

[14] G. T. Herman and A. Naparstek, Fast image reconstruction based on a Radon 
inversion formula appropriate for rapidly collected data, SIAM 3 . AppL Math., 
33, (1978) pp.511-533. 

[15] R. M. Lewitt. Reconstruction algorithinB : Transform methods, Proc. IEEE, 
71(3), (1983), pp. 390-408 

[16] P. Munshi, R. K. S. Rathore, K. S. Ram and M. S. Kalra, Error estimates for 
tomographic inversion, Inverse Problems,!, (1991) pp.399-408. 

[17] F. Natterer, The Mathematics of Computerized Tomography, John Wiley & 
Sons, New York, (1986) 


BIBLIOGRAPHY 


83 


[18] Kwoli, L S. Reed, T. K. Truong and C. M. Chang, 3-D reconstruction for diverg- 
ing X-ray beams, IEEE Trans. : Nucl. Sci., bf NS-25(3), (June 1978) pp.l(X)&- 
1014. 

[19] A. Madcovald, Physical problems of computerized tomography, Proceedings of 
the IEEE, 71(3), (1983) pp.373-378. 

[20] V. Bhatt, P. Munshi, and J. K. Bhattacharjee, Analysis of the performance 
of medical CT scanners using fractal theory. Bulletin of Radiation Protection, 
13(2), (1990) pp.38-40. 

[21] R- K. S. Rathoie, P. Munshi, and R. K. Jarwal, Measurement of void-fraction 
distribution using a tomographic chord segment inverrion algorithm, NucL Tech- 
noL, 82 (1987), pp.227-234 

[22] W. H. Miller. The design of a portable CAT scanner for utility pole inspection. 
Trans. American. Nuclear. Society, 52, (1986) pp. 350-351 



